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Abstract 

We exhibit a one to one correspondence between some universal probabilistic proper- 
ties of the ordering coordinate of one-dimensional Ising-like models and a class of contin- 
uous time random walks. This correspondence provides an new qualitative picture of the 
properties of the ordering coordinate of the Ising model. 

1 Introduction 

Many systems in Nature exhibit a phase transition when a control parameter (the reduced 
temperature for example) e is varied. Often, the phase transition is related to a spon- 
taneous symmetry breaking in the system: one phase is less symmetric than the other, 
and the so-called order parameter <p describes this symmetry breaking. The parameter <fi 
could be a scalar, a vector or a tensor and its number of components n is an important 
parameter. 

The amplitude of the order parameter is zero in the most symmetric phase and nonzero 
in the less symmetric one. If the order parameter varies continuously at the transition, 
one speaks of a second order phase transition. It turns out that, near a second order 
transition point, every physical quantity A (like the specific heat, the susceptibility, ...) 
can be written as the sum of a regular part and a singular one in e. The singular part 
(A s ) behaves as power law A s ~ \ e \~ a providing that the control parameter e has been 
chosen as to vanish at the critical point. The parameter a is called the "critical exponent" 
associated to A. 

A remarquable feature of second order phase transitions is the so-called "universal- 
ity". Scaling theory [1] and renormalization group arguments [2] show that, the critical 
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exponents are only function of n and d, the dimensionality of the system. Thus systems 
as different as for example an Ising model (describing the paramagnetic ferromagnetic 
transition in one component spin system) and a model of a fluid (describing the liquid gas 
transition) belong to the same universality class. 

A less investigated problem concerns the universal aspect of ordering coordinate un- 
derlying critical point singularities. A. Bruce [3] has investigated this problem in the 
case of an Ising-like problem, i.e. a problem for which the scalar order parameter is the 
ensemble average of a one-component local spin variable a(x). The configurational en- 
ergy consists of local double well potentials and the short-range interaction among the 
spins. Each stochastic variable a(x) is characterized by its probability density distribution 
P(cr). This function has a whole spectrum of possible forms, but two extreme cases can 
be singularized. At low temperature, the thermal energy is much smaller than the depth 
of the potential wells. Thus the ordering coordinate are localized in the vicinity of the 
two well minima and P(a) has a strongly two-peaked form. One speaks of the "ordered" 
limit. Alternatively, at high temperature, the potential wells are so shallow that the local 
coordinate is nearly distributed about the symmetric state a = 0. One speaks of the 
"disordered" limit. 

Parallel to the spin system, the continuous time random walks, (also called randomized 
random walks in W. Feller, [4]), do play an important role in the modeling of transport 
phenomena in disordered media, (see [5], for a recent bibliography on this topic). Firm 
mathematical basis to describe this class of processes, can be found in [6] and [7], the 
later being more directly related to one aspect of the models to be described here. The 
simplest situation consists of a symmetric random walks on the infinite one-dimensional, 
(1 — D), lattice with lattice spacing a and with the random time between consecutive 
the jumps being exponentially distributed. In this case, the position of the walker is 
described by a markovian Master Equation. When the time between the jumps is drawn 
from non-exponential distributions, non-markovian properties arise and to describe these 
situations E. Montroll and G. Weiss proposed a generalization of the Master equation [8]. 
In the present paper, we shall consider the non-markovian random walks obtained when 
the time intervals between jumps are drawn from generalized Erlang distributions. These 
distributions defined on [0, oo] belong to the general class of phase-type distributions which 
describe the time until absorption of a finite state continuous time Markov chain [9]. For 
this class of random processes, we shall observe that, when large time are considered, the 
probablistic properties of the position of the walker exhibit a universal behavior which is 
in one to one correspondence with the probabilistic properties of the ordering coordinates 
of a one-dimensional Ising spin systems. 

Our paper is organized as follows in section 2, a brief account of the universal properties 
of Ising-like systems are given. In particular, the probability distribution of a collective 
coordinates in an Ising chain is recalled. In section 3, we introduce the relevant class of 
continuous time random walks to be used and calculate the probability density describing 
the position of the walker. In the last section, we exhibit the one to one correspondence 
which can be drawn between the two classes of models. 

2 Universal properties of the probability density 
distribution of the collective coordinates in Ising- 
like models. 

Let us return to the Ising-like system described in the introduction. To be able to de- 
tect the universal features of the ordering coordinates, one has then to coarse-grained 



the problem and considers "bloc-spin" variables o\l(x) representing the effective spin of 
a block of linear size L. The variable &l is also a random variable characterized by its 
probability density distribution P(ctl, L). Renormalization group arguments suggest that 
Pl(^l, L) tends to an universal limiting form P* when both L and the correlation length 
£ are much larger than a, the lattice spacing. It turns out that in one dimension the prob- 
ability density function P* can be computed exactly using a matrix transfert operator 
technique [3]. Furthermore, the critical temperature of one-dimensional Ising-like systems 
with short range interactions is zero. The low temperature properties of the block coor- 
dinate probability distribution function, for large block size, is governed by the statistical 
mechanics of the cluster walls [10]. Thus to compute P{<jl,L) for our one-dimensional 
Ising like models, it suffices [3] to consider a one dimensional fixed-length spin Ising model, 
a(xi) = ±<7q, with nearest neighbors interactions and the following effective Hamiltonian 
(in units of ksT, where ks is the Boltzmann constant and T the temperature): 

K N 

£ = 2 <?(xi)<r(xi+i) (1) 

°o i= i 

<7o is chosen such as P(<tl, L) has an unit variance. 

It is known [11] that (in the thermodynamic limit N — > oo) the two-point correlation 
function of this model is: 

C 2 (r) =< a(x)a(x + r) >= al [tanh(K)] r = a\ exp(-r/£) (2) 

where 

i = l - ex P (2if )[1 + 0(exp(-4ff))] (3) 

is the correlation length. 

The block coordinate takes the simple form 

1 L 

°L = T— I>( X *)- ( 4 ) 
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To compute P{ai,L) it is suitable to introduce its Fourier transform P(H,L): 

1 f°° ~ 

P(a L ,L) = — P(H,L)exp(-iH<T L )dH, (5) 

^ J-oo 

P(H, L) is simply the characteristic function from which the n-point cumulants can be 
derived. Moreover, P(H,L) can also be written as [3]: 

W)^, (6) 
where Z(h) is the canonical partition function defined as 

h N 

Z(h) = Tr {a(Xi)} exp[-£ + — ]T a(x z )}. (7) 

The traces Z(h) and Z(0) can be readily computed by the usual transfert matrix method 
[11] for h real. The analytic continuation to the imaginary h axe is unique. 
In the large L and £ regime, one finds [3] that P(H, L) is indeed an universal function, 
denoted P*(H, z), depending only upon the reduced variable z = L/£: 




P*(H, z) = exp(-|) { cosh \\l ~— - H 2 \ + ; ^== sinh ( \\ ~— - H 2 \ } . (8) 
By inverse Fourier transform one finds: 




z , z 



P*(a L , z) = \ exp(-|){ (5(1 - a L ) + 5(1 + a L )) } + 

Ii(|(l-^) 1/2 ^ 



exp(~) I ( |(1 - alf 2 ) + 'V o :' /9 ' 6(1 - oj). (9) 



2 V L ' J d-alf 2 
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Two particular limits deserve a special attention: 

• i). The critical limit: z — > 0. 
In this case one finds: 

P*(a L , z)n^ l 5 ^ ~ °l) + 5(1 + <tl)] • (10) 

Hence the probability distribution function has a limiting "ordered" form. 

• ii). The non-critical limit: z — > oo. 
In this cas, one finds 

T r-zo\ 



which corresponds to the "disordered" limit. 

3 Semi- Markov Birth and Death processes 

Now we consider a symmetric birth and death process on the infinite 1 — D lattice: 
Z a = {..., —a, 0, +a, 2a, ...}, a being the lattice spacing. The jumps of the random walker 
occur at random times £o, £i> £2 and the interval of time between successive jumps is 
denoted = £k+i ~ £fc > k = 0, 1, 2, ... . We assume the to be independent and iden- 
tically distributed random variables. Let ip(\,jl,t) be the probability density from which 
the {t/J's are drawn. The constant vectors A = {Ao, Ai, ....Ajv} and fl = {/1q, ■■■Hm} 
are parameters characterizing the distributions of the time between consecutive jumps. 
Accordingly, we shall write: 

Prob {t < t < t + dt} = V>(A, fl, t)dt. (12) 

From now on, we shall assume that tp(\,(2,t) is called a PH-Type density, which can be 
defined by its Laplace transform: 







, n- 

r/.(A, fl, t)e~ ut dt = V(A, fl, u) = -f (13) 



with the definitions: 



n2U(l + £), foriV>l, 

n x = <; (14) 

1, for N = 0. 



and 



nf=i(l + ^), forM>l, 
% = \ (15) 
1, for M = 0. 

We impose that N > M. In view of Eqs.(|l3|), (|i~4l ) and (|l5|), we immediately obtain : 



and the variance of the time between consecutive jumps reads as: 



d 2 ~ - 

°"r = ( r2 ) - ( T ) 2 = Q^2^{\rf, U ) \u=0 



~-ip(X,fl,u) \ u= o 



N j ' 

Y — 



N 

E 



M -i "I M -. 

E: 1 +E- 1 



,fc=i 



(17) 



Let us now return to our random walk and denote by P(sa, t) the probability to find the 
walker at the position sa at time t. We shall assume: 



P(sa, 0) = 5 sat o and £ = 0. 



(18) 



where 5 sa fi is the Kronecker symbol. Using the formalism introduced in [5] and [6], it is 
shown in [12] that P(sa, t) obeys to the high order differential difference equation: 

N ^ f) M i f> 

n (i + TwJ p ( sa > *) = n (! + -757) t p + ^ *) + p (( s - ^ *)i • (19) 

fe=0 AfeC ^ fc=0 ^ fcCW 

To solve Eq.(19), we would need to specify the appropriate initial conditions, -j^P(sa, t) \t=o 



for m = 1, 2, iV— 1. As we shall not derive, in this paper, solutions of Eq.(|19|), we refrain 
to give here further details on the characterization of these initial conditions. 
Let us now give two simple illustrations: 



Example 1 : With N = 1 and M = 0, Eq.(|T|) takes the form: 



1 + 



J_ d_ 

X^dt 



P(sa, t) = -[P ({s + l)a, t) + P ((s - l)a, t)} 



(20) 



which is the Master equation for the continous time Markov random walk. 
Example 2 : With N = 2, A = Ai and M = 0, Eq.(^) takes the form: 

( d 2 d \ 



\2 
A 



-P(sa, t) + ~P ((* + l)a, t) + ~P ((a - l)a, t) 



with the initial condition being for this case: ^P(so, t) |t=o- 



(21) 



For small lattice spacing a, we expand the right hand side of Eq.(19) up to second order 
in a: 



N 



fe=l 



1 d 



M 



n( i +va)^>=n(i+ 



it=i 



P(x,t) + 



2 



:P( X ,t) 



(22) 



where we have written sa = x £ 1Z. As usual, we now rescale the time variable as: 

T = a 2 t. (23) 

Introducing Eq.(|2^) into Eq.(p2|), we obtain: 

N 



E C k a 2k 



P(x, T) + a 4 C 2 — ,P(x, T) + a 2 d— P(x, T) 



fe>2 



<9T fe 



<9T 



2 pp. M r, k o 2 

" FM)iE« a+2 4^.n ( 24 ) 



2 dx 2 v ' dT k dx 2 



where and Bk are constant coefficients Vfc. Using Eq. ([[(]), we have in particular: 

d = (r). (25) 



When M = 0, the PH-Type distribution Eq. (|13[) reduces to the E n -type distribution, 
(generalized Erlang distributions of order n). From now on, we shall focus our attention 
to this class of E n -type models. Retaining the terms up to order a 2 in Eq.(|2"4|), we have 
therefore when M = 0: 

a 2 C 2 ^- 2 P*(x,T) + (r)^P*(x,T) = \^^T). (26) 

where the superscript * for P*(x,T) indicates that we restrict ourselves to the E n -type 
model. In view of Eqs. (pi), fll6|) and (O), we can write: 



C2= 1 -((rf-al) = 1 -(r) 2 (l-CV T ^ (27) 

where CV 2 . stands for the coefficient of variation which for E n -type distribution is always 
smaller than unity. 

Eq,(|26|) is the celebrated Telegrapher's equation which, in the asymptotic regime, appears 
as a "universal" description for all random walks with time between the jumps drawn 
from E n -type distributions. 
Introducing the Fourier transform: 

P*(k,T)= P*{x,T)e- ikx dx, (28) 



the Telegrapher's Eq.(|26D takes the form: 



a 2 C 2 ^P*(k,T) + ( T )—P*(k,T) = -P*{k,T) (29) 



and the solution of Eq.(p9[) with initial conditions: 



P* ( k , T = 0) = and ^P*(k,T) \ T=0 = (30) 

reads immediately as in Eq.(^|), provided we reinterpret the parameters as: 

2T 



and H 2 = ^-z 2 k 2 , (31) 



q(r)a 2 4 
with 

q=(l-CY 2 T ) € [0,1]. (32) 
By direct inversion, the density Eq.(^) follows provided we introduce the identification: 



2 ^ /oq^ 

z z a z q 



Although this is not necessary, we can chose T = 1 which corresponds to the definite scale 
factor chosen in [3], (see Eq. (2.11) in A. Bruce, [3]). 

Hence, the critical and non critical regimes described by Eq,(|i~0|) and ( |Tl| ) follows directly. 



4 Conclusion 



The Telegrapher's equation Eq.(|26|) and its solution Eq.(^) describe the asymptotic be- 
havior of the probability density for both the coarse-grained variables of an Ising-like chain 
and the position of continuous time random walker. 

The behavior of the block coordinate of the Ising system is similar to the behavior 
of a random walker with time between the jumps drawn from an E n type distribution. 
This can be understood qualitatively as follows: consider the coarse-grained variable of 
the Ising model including a given, (large), number of spins, say iV >> 1. Add one spin to 
this sample of N spins and consider the new probability density distribution of the block 
spin coordinate. Two main factors influence the result: 

• 1) the correlation between nearest-neighbor spins. 

• 2) the fact that in the large sample N limit, only the spins with values close to ±o"o 
will contribute to the coarse-grained variables. 

Hence, to an increase of the sample used to construct the coarse-grained spin variable in 
the Ising model corresponds an increase of the time in the continuous time random walk 
with a lattice sapcing being related to aQ. Furthermore, the "aging" effect implied by 
the use of non-exponential distributions for the time intervals between the jumps directly 
describes the correlation between the spin variables in the Ising model. Quantitatively, 
the roles played by the parameters of the two models are as follows: to the parameter 
tanh(iT) 6 [0,1], K > defined in Eq.(^), we can directly associate the parameter 
q G [0,1] defined in Eq.(|32|). This correspondence is physically consistent as tanh(iT) 
measures the correlation between nearest-neighbor spin variables given by Eq.(|2|) while 
q relates the memory effects of the semi-markovian random walk. Two limiting regimes 
deserve attention: 

• a) Weak correlation regime. Here q — > <^ CV? = 1 and hence tanh(X) 44> 

K — > 0. We then reach here the markovian regime of the random walk. Accordingly, 
the z — > oo limit arises in Eq.(|3l|) which implies an absence of correlations between 
the spins. 

• b) Strong correlation regime. To the limit K — > oo tanh(i^) — > 1 corresponds 
the q = 1 => CV^ = regime for the continuous time random walk. As CV^ = 
0, the time intervals between the jumps become equally spaced and this therefore 
corresponds to a discrete, (deterministic), time random walk. Accordingly, we have 
z — > and therefore the double-peak shape of the probability density exhibited in 
Eq.(|l0|) arises consistently. 
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